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Abstract 

Many of graphene's unique electronic properties emerge from its Dirac-like electronic energy 
spectrum [IHZ!- Similarly, it is expected that a nanophotonic system featuring Dirac dispersion 
will open a path to a number of important research avenues [9l-ll4j. To date, however, all proposed 
realizations of a photonic analog of graphene lack fully omnidirectional out-of-plane light confine- 
ment, which has prevented creating truly realistic implementations of this class of systems. Here 
we report on a novel route to achieve all-dielectric three-dimensional photonic materials featuring 
Dirac-like dispersion in a quasi-two-dimensional system. We further discuss how this finding could 
enable a dramatic enhancement of the spontaneous emission coupling efficiency (the /3-factor) over 
large areas, defying the common wisdom that the /?- factor degrades rapidly as the size of the system 
increases. These results might enable general new classes of large- area ultralow-threshold lasers, 
single-photon sources, quantum information processing devices and energy harvesting systems. 
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Since its isolation from bulk graphite in 2004 [T], graphene -a one-atom-thick sheet of 
carbon- has attracted an ever increasing amount of interest [2HZ]; nowadays, the study 
of the electronic properties of this two-dimensional (2D) material has become one of the 
most active areas of condensed matter physics. This general endeavor has also stimulated 
new directions in related research fields, especially those originally inspired by the physics 
of electronic transport in crystalline solids. Of particular relevance in this context are 
photonic materials whose dielectric constant is periodically structured at the subwavelength 
scale, the so-called photonic crystals (PhCs) [8J. By exploiting the analogy between the 
propagation of electrons in graphene and the propagation of photons in suitably designed 2D 
PhCs, phenomena such as directional optical waveguiding [9], pseudodiffusive transport of 
light p^ITT]. Klein tunneling [12], and the observation of the Zitterbewegung of photons [13] 
have been recently proposed. The existence of Dirac points at the center of the Brillouin 
zone induced by accidental degeneracy in square lattice 2D PhCs, has also been discussed 
recently [13]. However, all these systems share the common fundamental drawback that 
they lack fully omnidirectional out-of-plane light confinement, which has so far prevented 
the creation of a truly realistic implementation of a photonic counterpart of graphene. 

In this Letter, we propose a feasible approach to achieve simultaneously quasi-two- 
dimensional light propagation and Dirac cone dispersion in an all-dielectric 3D photonic 
material particularly suitable for optical device integration. We show how the unique light 
confining properties of a proper choice of 3D layered PhCs enables creating extended planar 
defect modes whose dispersion relation exhibits isolated Dirac points inside a complete 3D 
photonic band-gap. In the limit in which the emitter frequency virtually coincides with the 
Dirac point frequency (i.e., the frequency of the Dirac cone vertex) the number of photonic 
states available to the emitter approaches one, even if the system features a macroscopic 
area. Thereby, the photonic materials presented in this Letter enable for the first time the 
implementation of structures much larger than the wavelength, which nevertheless have f3 
factors close to unity. Due to the crucial role played by the /3-factor in various areas of 
physics (from optoelectronics to quantum computation or energy harvesting), we believe 
that these results hold a great promise for the development of novel types of nanodevices. 

A schematics of the considered photonic material is rendered in Fig. la. We start by 
considering the electromagnetic properties of the PhC before the defect plane is introduced; 
it is a face-centered cubic (fee) 3D PhC of air (or low-index) cylinders embedded in a 
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FIG. 1. An isolated photonic Dirac point in a 3D photonic crystal. Panel (a) displays a 
sketch of the considered structure. The basic structure consists of face-centered cubic (fee) lattice of 
overlapping air cylinders embedded in a high dielectric background (yellow regions). This structure 
can also be seen as an alternating stack of rod and hole layers (see labels in the figure). A planar 
defect is introduced in the structure by removing a hole layer and replacing it by a triangular array 
of finite- height dielectric cylinders (green cylinders). The geometrical parameters used to define 
the structure are also displayed. The inset renders the electric-field intensity corresponding to a 
guided mode in the defect plane at the Dirac frequency. Panel (b) displays the dispersion diagram 
corresponding to the structure shown in Panel (a), projected over the First Brillouin zone of the 
in-plane triangular lattice characterizing the hole and rod layers of that structure. The inset of 
Panel (b) shows an enlarged view of the dispersion diagram near the Dirac point. 

dielectric background and oriented along the (111) direction [T5| 116]. This peculiar class of 
layered PhCs can actually be viewed as an alternating stack of two different types of layers. 
One of the layers has the form of a triangular lattice of finite-height dielectric rods in air 
(labeled as rod layer in Fig. la), whereas the other layer can be described as a triangular 
lattice of air holes milled in a dielectric slab (labeled as hole layer in Fig. la). This PhC 
features two important characteristics. First, each of the two types of layers displays a 
highly symmetric cross section that mimics a canonical 2D photonic-crystal structure: one 
is a periodic array of air holes in a dielectric dielectric slab, and the other is a periodic array 
of hexagonal-like rods in air. Second, although neither of the layers displays a complete 
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(omnidirectional) photonic band gap by itself, when the layers are periodically stacked as 
shown in Fig. la, a large complete photonic band-gap can be obtained, using practical values 
of the refractive index contrast [TSl US]. 

Next, consider creating an extended planar electromagnetic defect mode in the 3D PhC 
described above. In order to that, we introduce into the structure a single defect layer 
that perturbs the original periodic sequence rod-layer/hole-layer/rod-layer along the (1 1 1) 
direction. In particular, we remove a hole-layer of the structure and replace it by a triangular 
lattice of finite-height dielectric rods with circular cross section (green cylinders in Fig. la), 
whose radius, height and dielectric constant are given by r^, hd and e^, respectively. 

It is known |16l [T7] that the introduction of line defects into layered 3D PhCs of the 
type described above enables the implementation of localized electromagnetic states whose 
dispersion relation, field profiles, and polarization are in a close correspondence with those 
associated with the corresponding 2D PhC geometries. Similarly, one would expect that 
the spectral properties of the planar extended states localized in the defect layer rendered 
in Fig. la should, to some extent, inherit properties associated with the Bloch states of the 
bona-fide 2D counterpart of this defect layer (i.e., those that are present in a triangular array 
of infinitely long high-index rods). On the other hand, it has been demonstrated that the 
intrinsic symmetry properties of 2D PhCs based on a triangular lattice can induce the pres- 
ence of Dirac points near high-symmetry points of the band structure [914TT] . Thus, this line 
of reasoning suggests the feasibility of creating a single 3D physical system simultaneously 
featuring quasi-two-dimensional light propagation and Dirac cone dispersion. In addition 
to these two features, it is also crucial that the mentioned Dirac cone dispersion is isolated 
within a given frequency bandwidth, or, equivalently, that the Dirac cone is fully separated 
from the rest of the bands present in the band structure. Only by combining all three of these 
features in the same system, it is possible to fully exploit the analogy between electronic 
and photonic graphene. We emphasize that, although, in principle, it is not challenging to 
obtain each of these three characteristics separately in a photonic structure, achieving all 
three of them concurrently in the same 3D physical system is not straightforward. 

To explore to which extent these ideas can be implemented in a realistic system, we 
have carried out a systematic numerical analysis of the evolution of the corresponding band 
structures as a function of the geometrical parameters and the dielectric constants, defining 
both the underlying 3D PhC, as well the defect layer. The calculations were performed 
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by means of the plane-wave expansion method to Maxwell's equations [18] using a supercell 
large enough in the (1 1 1) direction so the properties of an isolated defect plane in an infinite 
3D PhC are accurately reproduced. Our calculations show that, for the optimal structure, 
the air holes within the hole layer and the equivalent- cylinders [IT] in the rod layer are 
Vh = 0.41a and Vc = 0.18a respectively (a is the lattice constant of the in-plane triangular 
lattice defined within each layer; see Fig. la). The thicknesses of the hole and rod layers are 
th = 0.32a and = 0.50a, respectively, whereas the refractive index of the high-dielectric 
material is assumed to be n = 2.5. The low- refractive index of the structure is taken to be 
air. On the other hand, the defect layer of this optimal system features = 3.1, = 0.32a, 
and hd = t^. In addition, from this numerical study we have also obtained that the optimal 
configuration is that in which the cylindrical rods of the defect layer are aligned with the 
rods forming the two rod- layers located right above and below the defect plane (see Fig. la). 

Figure lb summarizes the dispersion diagram obtained for this structure. Shaded violet 
areas in this figure show the projected band structure for the perfectly periodic 3D PhC 
(i.e., without the defect layer). In this case, the dispersion diagram was obtained by plotting 
the frequencies u of the extended bulk states of the system as a function of the in-plane 
wavevector ky in the irreducible Brillouin zone of the underlying in-plane 2D triangular 
lattice. The considered system exhibits a large 3D complete photonic band gap (shaded 
yellow area), centered at frequency uj = 0.497(27rc/a) (c is the light velocity in vacuum) 
and featuring a gap-mi dgap ratio of approximately 8%. The dispersion relations of the 
guided modes of the defect plane are also rendered in Fig. lb (red lines). As can be seen, 
these defect bands indeed display a Dirac point at ud = 0.506(27rc/a) (see Fig. lb and the 
inset of Fig. la), fully lying within the omnidirectional photonic-band gap of the periodic 
system (yellow area in Fig. lb). Importantly, this Dirac point is completely isolated from 
all of the rest of the frequencies of the band structure of the system within a bandwith 
Au = 0.026(27rc/a). Finally, the strong out-of-plane photonic-band gap confinement of the 
electromagnetic fields at u = ud is clearly observed in the inset of Fig. la, which displays 
the corresponding cross-section along the xy plane of the electric-field intensity distribution. 
We believe that these results represent the first time that a complete photonic analog of 
graphene is implemented in a realistic 3D physical system. We emphasize that the obtained 
results are scalable to many different frequency regimes, and therefore, could be used to 
enhance performance of different classes of active optical devices. 

5 



We now turn to demonstrate how the proposed class of photonic systems can enable an 
unprecedented control of light-matter interaction over large areas. In order to do that, we 
apply the commonly employed procedure to probe the properties of light-matter interac- 
tion in a complex electromagnetic environment: we study the radiation process of a point 
quantum emitter embedded in the system under analysis [3], [TT] . The emitter is modeled 
as a two-level system, characterized by a transition frequency Ug, an emission bandwidth 
(i.e., a transition bandwidth) of Aug, and a dipolar transition moment d = dd. In partic- 
ular, we analyze the spontaneous emission coupling efficiency, the so-called /3-factor. In a 
given multimode photonic system, the /3-factor quantifies the portion of all spontaneously 
emitted photons that couple into a certain targeted mode [22]. This physical magnitude is 
of a paramount importance in modern optoelectronics and quantum information processing 
since the performance of an active nanophotonic device, or a single photon emitter, can 
often be greatly enhanced by increasing the value of /3 [22H23j (recent examples include 
ultralow-threshold lasers [23H2Z] and single photon sources |2H1 EH] based on PhC cavities). 

The dependence of the /5-factor on the particular electromagnetic environment in which 
the considered emitter is embedded can be elucidated by examining its link with the corre- 
sponding photonic local density of states (LDOS), p{r,d,ijj). For a non-dissipative system, 
the LDOS can be written as [3] 

p(r, d,Lo) = J2 e(r) |E,(r) ■ dp 5{lo - u,) (1) 

where the index v labels the different source-free normal solutions to Maxwell's equations 
obtained for the considered photonic structure; Ei,(r) and lOy are their corresponding i?-field 
profile and frequency, respectively. e(r) stands for the dielectric constant distribution. 

On the other hand, from Fermi's golden rule [31 HU |19], one finds that, in 3D, the 
spontaneous emission rate of the considered quantum emitter F is proportional to the LDOS 
accessible to the emitter F = (Trldpcj^/S/ieo) pijs-, d, Ws). Thus, assuming the targeted mode 
to be a normal mode of the system, Et(r), of frequency w^, one can write the /3 factor as 
(see Suppl. Information), 

^ ^ ujtg{oJt)e{Y) |Ef(r) - dp 
/rfw5((a;) a;p(rs,d,w) 

where giiS) is the lineshape of the transition, centered at bjg and characterized by full-width- 
half-maximum (FWHM) of Acj^. Note that the factor giiot) in the numerator accounts for 
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FIG. 2. Dispersion relations of three representative photonic materials and typical 
spectrum of a quantum emitter. Panel (a) and (b) sketch, respectively, the dispersion relations 
of a homogeneous material and a photonic-crystal exhibiting a photonic band gap near the peak 
of the emission cOg- Panel (c): same as (b) but for the case in which the structure exhibits a Dirac 
point near ajg- Panel (d) renders the lineshape of the emission spectrum g{uj) in the form of a 
Lorentzian centered at cOs and featuring a full-width at half- maximum (FWHM) of Aug. Ug and 
ujd denote the central frequency of the band-gap and the Dirac point frequency, respectively. 

the fact that (3 decreases as the emission frequency is detuned from cot- Hereafter Ut = ujg is 
assumed. 

Equation (2) clearly shows that the /3 factor can be enhanced by introducing a phys- 
ical mechanism that minimizes the density of photonic states lying within the transition 
linewidth. In fact, the large values for /3 reached in subwavelength volume photonic res- 
onators [25], EHl EHl [29] can be viewed as a particular instance of this physical picture. Such 
nanoresonators are designed to have a volume small enough so that only one resonant mode 
lies within the transition linewidth. This makes them to effectively act as single-mode struc- 
tures, which as deduced from Eq. (1) leads to values (3^1 (provided that the coupling with 
the radiation modes existing outside the resonator is negligible). Although very relevant in a 
number of contexts, this cavity-based approach does not admit a straightforward extension 
for large-area control of (3. In the rest of this Letter, we show how the general class of 
isolated Dirac points discussed above do enable such large-area control. 

To gain physical insight into the effect of the Dirac cone dispersion on the (3 factor, we 
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evaluate the magnitude of ^ in the following three cases (all of them 3D): first, the case of a 
homogeneous material; second, the case in which the dielectric material is periodically struc- 
tured so that 0;^ lies in the vicinity of the lower edge of a 3D photonic-band gap; and, third, 
the case of the structure displayed in Fig. 1, i.e., a system exhibiting simultaneously quasi- 
two dimensional light propagation and an isolated Dirac point near Ug- Figures 2a, 2b, 
and 2c illustrate each of these dispersion relations (for simplicity in the visualization, in 
Figs. 2a and 2b, only the 2D counterparts of the corresponding cases are rendered). Impor- 
tantly, in all three cases we also assume that the EM field in the system is confined in a 
finite volume V, such that V » \^ (the dependence on volume is addressed below). 

The homogeneous case is characterized by the following dispersion relation: u:{k) = 
ck/n (where n is the refractive index and k = |k|, with k = (k^jk^))- Then, making 
use of the isotropy of the medium, one can derive the following simple expression for the 
P factor: ^ = (l/V) uJsg{uJs)/Fh{uj), where Fh{uj) = j du g{u) u ph{u), with ph{u) = 
(l/27r)(n/c)^a;^. Similarly, for the case of the photonic band-gap, taking ijj{k) — Ug — Agk'^ 
{ug being the center of the gap and, see Fig. lb, and Ag is a constant which has to be 
determined from calculations of the band structure; physically, Ag defines the curvature of 
the dispersion relation close to Ug), one finds that the magnitude of /3 can be calculated 
using the same expression given above for the homogeneous case but replacing Fh{uj) by 
Fg{uj) = Jduj g{uj) CO pg{uj), with pg{uj) = {A7r^Af^)-\ujg - a;) V^. 

For the Dirac case, the calculation of (3 is more involved that in the previous two cases. 
First, we take into account that an excited dipole embedded in the defect layer displayed in 
Fig. la only decays via the guided modes confined within the layer; the decay into any other 
modes surrounding the layer (e.g., bulk Bloch modes) is prevented by an omnidirectional 
photonic band gap. Then, for small enough values of ha (so only the fundamental mode 
guided mode in the 2;-direction is excited), the quasi-2D light propagation inside the defect 
layer can be described by the dispersion relation corresponding to in- plane Bloch states (i.e., 
states with = 0): w(k||) = Ud ±A£,\k\\ — ky o| (see details in Suppl. Information). In this 
expression, A^ is a constant that can be obtained from band-structure calculations and k||_o 
defines the coordinates in k-space of the Dirac cone vertex, whereas the plus and minus signs 
correspond to a; > ujd and ou < ujd, respectively. (Note that, physically. Ad corresponds 
to the slope of the Dirac cone). Thus, one finds that the total density of states accessible 
to the emitter as Pd{'^) — l/(27rA|,)|a; — ud\- This, in turn, yields the following expression 
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for the SE coupling efficiency ^ — {1/ A)ujsg{ujs)/ Fd{oS)^ witfi Fd{uj) — J dujg{uj)u;pD{co). 
Here A is the transversal area of the defect layer (i.e., the total volume of the defect layer is 
V = Axha). 

Wc now quantify the values of (3 for each of the above cases using realistic parameter 
values. To allow for a direct comparison among the three considered class of systems, we 
introduce a renormalized spontaneous emission coupling efficiency 77. In the homogeneous 
and band-edge cases we define this magnitude as r)h,g — Ph,g x V/a^, whereas for the Dirac 
case we define rjo — x S/a^ . This normalization allows us, on one hand, to remove 
from the discussion the obvious dependence of 13 on the size of the system, and thus, focus 
exclusively on the photonic properties. On the other hand, it also removes the geometrical 
factor V/A that enhances the /3 factor in the Dirac case which respect to the other two cases. 
This factor stems from the electromagnetic confinement in the direction of the guided 
modes in the Dirac structure, and therefore cannot be ascribed to the Dirac spectrum. 
Furthermore, in our calculations we assume that the transition lineshape is described by a 
Lorenztian centered at cu^ = 2.1 x 10^^ s~^ (i.e., an emission wavelength of 900 nm) and 
featuring a relative FWHM Aug/oOs = 10^^; values more than one magnitude smaller for 
AcUs/cUg can be reached using, for instance, quantum dots at low temperatures. In the 
homogeneous case a refractive index n — Udis chosen, whereas for the band-edge and Dirac 
cases, the values of the dispersion relation parameters Ag and A^ are taken from band 
structure calculations: Ag — 1.2 x ca/(27r) and Ad = 0.3 x 1/c (a is lattice constant defined 
in Fig. la, which for operation at the considered emission wavelength, takes a value of 450 
nm). 

Thus, using these parameters, from the numerical evaluation of the expression for (3 given 
above, we obtain rjh — 68.2, rjg — 4.0 x 10^ and rjo — 3.6 x 10^. As readily deduced from 
these values, the Dirac dispersion introduces an enhancement factor of about four orders of 
magnitude with respect to the homogeneous case and about one order of magnitude with 
respect to the band-edge case. This is an important result, since, as wc show below, enables 
reaching values of /3 ~ 1 over macroscopic areas. Physically, the origin of this dramatic 
increase of rj (and consequently of /3) can be understood in terms of the rapid decrease of 
the number of photonic states available to the emitter as its emission frequency approaches 
the frequency of the Dirac point. In particular, in contrast to the homogeneous and band- 
edge cases, in the Dirac case when oug — > ujd (i.e., as emission frequency approaches the 
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FIG. 3. Finite size effects on the /3-factor. Dependence of the /3-factor on the normahzed 
emission linewidth Aojg/uJs, as computed for several sizes of the transverse area A of the system, 
with a being the periodicity of the in-plane triangular photonic crystal (see definition in Fig. la). 

Dirac vertex frequency), the number of modes accessible to the emitter approaches unity, 
making the whole structure to effectively behave as a single-mode system, even if it features 
a large area. Note that the LDOS is strictly zero at Ug and in the band-edge and Dirac 
cases, respectively. Therefore, in the calculations for each case, we have assumed Ug to be 
slightly detuned (by a quantity much smaller than cOs and Aug) from Ug and ud- 

We now focus on analyzing the dependence on size of the enhancement of /3 enabled 
by isolated photonic Dirac cones. The physical origin of this dependence stems from the 
fact that for values of the area A such that A » cannot be safely assumed, it is 
necessary to account for the discreteness of the eigenmodes in the system: the only states 
that will be allowed to exist are those characterized by a wavevector k = {kx,ky) whose 
value coincides with one of the nodes of the rectangular grid defined by the discrete set of 
values {27mx/L,27Tny/L} (with and ny being arbitrary integers, and where the system is 
assumed to be square shaped with side length L, i.e., A = L^) [6]. Therefore, the influence 
of these finite-size effects on (3 can be computed by using the discrete version of Eq. (2) (see 
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Suppl. Information). 

Figure 3 shows the computed results for /3 as a function of the normahzed emission 
bandwidth Aug/us for several values of the lateral size of the system, ranging from L/a = 10 
to L/a = 10^. As observed, in the considered cases, /3 tends to 1 when Aus/oJs approaches 
the lower limit of the interval displayed in the figure {Auog/us = 10^^). This is due to the 
fact that for all the considered system sizes, a linewidth Aug/uJs ^ 10^^ is smaller than the 
frequency interval between the adjacent modes, and therefore, the structure is actually acting 
as a single mode system (much in the same way as occurs in large-/? photonic nanocavities). 
As Aus/us is increased, a growing number of modes enter into the interval where g{u) is 
not negligible and, therefore, the value of /3 starts decreasing. Since the frequency interval 
between adjacent modes is smaller for larger values of A, the decrease of the /3 factor with 
AlOs/oOs starts sooner for larger values of A. 

In conclusion, we have demonstrated for the first time that a 3D photonic material can 
exhibit isolated Dirac cones in its dispersion relation. The proposed approach is readily 
accessible experimentally, both for fabrication and measurement. In addition, we have shown 
how the proposed class of systems enables achieving unprecedentedly large values of /3-factor 
over large areas. We believe that these results could open new avenues in a variety of different 
fields. Thus, the large /3 factors demonstrated in this work may lead to the realization 
of important applications such as low-threshold large-area lasers, enhanced single-photon 
sources and novel efficient platforms for solar energy harvesting. 
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Supporting Information for: Enabling single-mode behavior over 
large areas with photonic Dirac cones. 



In the following, we present a detailed description of the theoretical treatment of the 
spontaneous emission coupling efficiency (the P factor) used in the main text. We start by 
deriving a general expression for the /3 factor for the class of structures considered in our 
work. Within the Wigner-Weisskopf approximation [T], the rate of spontaneous emission F 
of a two-level quantum emitter embedded in a complex electromagnetic (EM) environment 
can be described by Fermi's golden rule [21 [3] 

r = ^^y^p(rs,d,a;,) (1) 

where d = c? d is the dipolar moment associated to the radiative transition of the emitter. 
Here Ug is the frequency of that transition, whereas eo and h are the vacuum permittivity 
and the reduced Planck's constant, respectively. The scalar function p(r,d,c<;) represents 
the photonic local density of states (LDOS) accessible to the emitter at position r = r^. 

On the other hand, in the case of a non-dissipative system, in which the electric field 
can be expanded in terms of a complete basis of transverse orthonormal modes {Em(r)} of 
frequencies {cOm}, the LDOS can be expressed as [3H5] 

p(r, d,uj)=Y^ 6{u - e(r) |d E„(r) ^ (2) 

m 

where e(r) is the position dependent dielectric constant characterizing the system. Each of 
the modes in Eq. (2) satisfies the following orthonormality condition 

J dre(r)E„(r)E:(r) = 5^„ (3) 

with 6nm standing for the Kronecker's delta. The transversality condition reads 

V ■ [e(r) E^(r)] = (4) 

Note also that each mode profile Em(r) can be obtained by solving the following wave 
equation 

V X [V X E„(r)] = poe(r)u;^E„(r) (5) 
where fiQ is the vacuum permeability. 
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Now, by definition, the /5- factor can be calculated as /5 = Tt/Taii, where Tt is the spon- 
taneous emission rate into a given targeted mode (often a laser mode) and Taii is the total 
spontaneous emission rate into all the modes of the system (including the targeted one) [7] . 
Thus, assuming that the lineshape of the emission transition is defined by function g{uj), by 
inserting Eq. (2) into Eq. (1) and integrating the resulting expression over u, we obtain the 
following expression for Taii, corresponding to an emitter located at r = 

^aii = ^g^l^''^ J du g{uj) u f{uj, Vs) (6) 

where function f{u, r) is defined as 

/(a;,r) =^5(a;-a;™) |d-E„(r)|2 (7) 

m 

Equations (6) and (7) summarize well the physical origin of the total spontaneous emission 
decay in the considered system: on the one hand, the different terms in the summand 
of Eq. (7) account for the different modes in which a single frequency component u of 
considered emission transition can decay to. On the other hand, the integral in u appearing 
in Eq. (6) accounts for the continuous sum of these possible radiative decay paths for all 
frequency components of the emission transition. Note that, as expected, the lineshape of 
the emission, g{u), acts as a frequency dependent weight in this sum. The additional factor 
u multiplying g{uj) in the integral of Eq. (6) comes just from the proportionality factor that 
links the spontaneous emission rate and the LDOS (see Eq. (1)). 

This physical picture of the decay process also allows us to obtain an expression for 
Tt simply by singling out the contribution to Taii that stems from the considered targeted 
mode. In particular, if we define Ef(r) and Ut to be the targeted electric-field profile and its 
corresponding frequency, respectively, the magnitude of Tt can be obtained by substituting 
g{uj) by gt{uj) = 5{uj - Ut) in Eq. (7) and by replacing /(a;,rj by ft{uj,r) = 5{uj - uJt)\d ■ 
Ej(r)|2. This yields 

Tt = ^^^^ootg{u:t)\d.Et{r,)\' (8) 

Note that dividing Eq. (6) by Eq. (8), and using the definition of the LDOS given in Eq. (2), 
we recover Eq. (1) of the main text. 

We now focus on the application of the above formalism to calculate the /3 factor in the 
case in which the emitter is embedded in a three-dimensional photonic crystal (PhC). We 
assume that the considered PhC is characterized by a finite volume V = L^x LyX (where 
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Ly, Lz are the dimensions of the PhC along the x, y and z axis, respectively). Here 
we emphasize that in our theoretical calculations, this three-dimensional analysis is applied 
to the homogeneous and band-edge cases discussed in the main text. (The homogeneous 
case can be trivially considered as a periodic system with an arbitray periodicity). For the 
Dirac case, however, due to the out-of-plane subwavelength confinement of the EM fields 
introduced by the full photonic band gap, the analysis is performed in terms of the in-plane 
transverse area of the system, A = x Ly (see discussion in the main text). 

To analyze the finite-size effects on the /3-factor, without loss of generality, we assume the 
volume V (or transversal area A for the Dirac case) to be surrounded by Born-von-Karman 
boundary conditions (i.e., periodic boundary conditions our theoretical analysis admits a 
straightforward generalization to other types of boundary conditions). In this case, the index 
m used above to label the modes can be identified with {n, k, cr}, where n is the band-index, 
k is the wavevector of each Bloch mode (k lies inside the First Brillouin Zone (FBZ)) and a 
labels the polarization (cr=l and cr=2, for s- and p-polarization, respectively). In addition, 
since the system is finite, k can only take discrete values: k = 27r x [n^/ L^^Uy/ Ly^Uz/ L^) 
for the homogeneous and band-edge cases, and k = 27r x {rix/Lx., ny/ Ly^ 0) for the Dirac case 
(in all three cases, n^, ny, and are arbitrary integers). Thus, once the normal modes of 
the system E„^k,(7 are computed (for which we have used the plane-wave expansion method 
to Maxwells equations[8]), from Eqs. (6) and (7) the /3 factor can be calculated from 



In the limit in which the volume V of the system is such that >> A'^ (or equivalently, 
for the Dirac case, when the area A is such that A » A^), with A being the central emis- 
sion wavelength, semi-analytical expressions for the /3-factor can be obtained by assuming 
a continuous distribution of k vectors over the FBZ. Specifically, we can replace X^k ^ 
VI{2t{Y jp^^dk in the homogeneous and band-edge cases, and X^k ~^ ^/(27r)^ /p^^ (ik in 
the Dirac case (note that in all three cases the integral over k is performed over the whole 
FBZ) . Then, if we expand the argument of the Dirac delta appearing in the denominator 
of the right-hand side of Eq. (9) using 



and neglect the contribution of second order terms in |ko — ko(c(;„_o-)| [6], after some straight- 



ujt g{ujt) |Et(r) ■ d 





uj - u;„,a(ko) = Vkot^[ko - ko(wn,a)] + Odko - ko(u;„,<^)n 



(10) 
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forward algebra, one finds that Eq. (9) can be rewritten as 



1 ojsgj^^s) 

V J du giu) uj p{uj) 



(11) 



where the function p{uj) determines the total density of photonic states per unit volume in 
the structure. Note that in the Dirac case, V must be replaced by the transversal area A. In 
analogy with standard analyses in solid-state physics [U] , in the homogeneous and band-edge 
cases, p{co) can be expressed as 



where A{us) denotes the equifrequency surface u = cOs, and Vg is the magnitude of the 
group velocity Vg = \dco/dk\. At each point of the equifrequency surface, kt stands for the 
component of the 3D vector k that lies along the tangential direction to A{u)s) at each point 
of the fc-space. In the Dirac case a similar expression holds for p{uj), but now the domain 
of integration in Eq. (12) is a equifrequency curve instead of equifrequency surface. The 
resulting expressions p{u), obtained by performing the integral defined by Eq. (12) for the 
different dispersion relations considered in this work, are discussed in detail in the main 
text. 

Importantly, in deriving Eq. (11) we have assumed that for the range of parameters con- 
sidered in this work, the emission bandwidth is narrow enough so that |E„ko-(i")P ~ |Et(r)p 
for all modes whose equifrequencies Unka lie inside the interval where g{u) is not negligible. 
To verify numerically the accuracy of this approximation in the Dirac case (a similar analysis 
holds for the band-edge case), we have probed directly the LDOS of the 2D counterpart of 
the defect layer structure shown in Fig. la of the main text. Specifically, in order to do 
that, we have computed the power radiated by a dipole placed in the low-refractive index 
regions of the structure (i.e., in the interstitial regions among cylinders), Pout(i', d, w). To 
compute Pot(t(r, d, w) we have employed a generalization of the conventional coupled-mode 
theory [U] , in which each Bloch mode is considered as an independent input / output channel 
(see details in Ref. [IDj)- The computed results are summarized in Fig. 1 of this Supple- 
mentary Information, in which band structure calculations (Fig. la) are displayed together 
with the corresponding dependence of Po«t(i', d, w) with frequency (Fig. lb). As observed 
in these results, the dependence of Po„f(r, d, w) with co near the Dirac frequency ud (and 
hence the LDOS [H]), obtained by assuming |E„ko-(r)p ^ |E<,(r)p, is in good agreement 




(12) 
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Normalized Emitted Power P(r^,co,uJ 

FIG. 1. (a) Numerical calculation of two-dimensional photonic bands displaying an isolated Dirac 
point. The analyzed system is formed by a triangular lattice of dielectric cylinders of refractive 
index = 3.1 and radius = 0.32a (a being the lattice constant) embedded in air. (b) Power 
emitted by a dipole located at the center of the unit cell of the considered photonic crystal and 
with its dipolar moment pointing along the j;-axis (see sketch in top inset of this panel). In both 
panels, ud marks the frequency of the Dirac point. Bottom inset in Fig. lb shows the comparison 
between the predictions of semi-analytical and full numerical calculations (cyan and blue lines, 
respectively) for the emitted power. 



with full numerical calculations within a moderately large bandwidth of frequencies (see 
comparison between cyan line and blue line in the bottom inset of Fig. lb). Finally, we note 
that although in the particular case considered in Fig. lb the dipole has been placed at the 
center of the unit cell of the triangular lattice, with d pointing along the x-direction (see 
top inset of Fig. lb), we have checked that similar good agreement between numerical and 
and semi-analytical results is obtained for other positions of the dipole in the unit cell, as 
well as for other orientations of d. 
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